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ABSTRACT 



The time-dependence of heating in solar active regions can be studied by 
analyzing the slope of the emission measure distribution cool-ward of the peak. 
In a previous study we showed that low-frequency heating can account for 0% 
to 77% of active region core emission measures. We now turn our attention to 
heating by a finite succession of impulsive events for which the timescale between 
events on a single magnetic strand is shorter than the cooling timescale. We refer 
to this scenario as a "nanoflare train" and explore a parameter space of heating 
and coronal loop properties with a hydrodynamic model. Our conclusions are: (1) 
nanoflare trains are consistent with 86% to 100% of observed active region cores 
when uncertainties in the atomic data are properly accounted for; (2) steeper 
slopes are found for larger values of the ratio of the train duration A H to the 
post-train cooling and draining timescale Ac, where Ah depends on the number 
of heating events, the event duration and the time interval between successive 
events ijc); (3) tq may be diagnosed from the width of the hot component of 
the emission measure provided that the temperature bins are much smaller than 
0.1 dex; (4) the slope of the emission measure alone is not sufficient to provide 
information about any timescale associated with heating - the length and density 
of the heated structure must be measured for Ah to be uniquely extracted from 
the ratio A H /A C - 



Subject headings: Sun: corona 
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Introduction 



It remains a difficult question whether the solar corona, at temperatures exceeding 
a million degrees, is heated steadily or impulsively. Any potential heating mechanism 
must explain observed extreme ultraviolet (EUV) and X-ray emission, as well as observed 
dynamical activity. We must understand the mechanism, how it stores and releases energy, 
how that energy release affects the plasma and pr oduces the emissi on, and finally we must 



Klimchuk 



20061 ). Although we cannot 



predict observable properties of coronal loops (see 
directly observe coronal heating, we can study the time dependence of the heating which 
may give insight into the mechanism by which energy is released. 

Impulsive heating is character ized by short bursts of ene rgy release (which we refer 



to as nanoflares; see, for example, 



Cargill fc Klimchuk 



19971 ) and a period of cooling 



between successive heating events. If the cooling period is short then the plasma will not 
cool much between heating events and the heating may essentially be treated as steady 
(the limit of high-frequency nanoflares), resulting in a more or less isothermal plasma 
distribution. Contrariwise, if successive bursts are sufficiently far apart then the plasma 
will cool significantly and, at any particular time, plasma on different magnetic strands of a 
multi-stranded loop may have a broad distribution of temperatures. Thus, the temperature 
distribution of plasma, quantified by the emission measure (EM), may help to distinguish 
between different heating scenarios. We will focus on heating in active region cores in the 
present work. 



As shown in the first paper of this series 



Bradshaw et al. 



(120121 ) (hereafter referred 



to as Paper I), EMs in active region (AR) cores can be characterized by the slope a (i.e. 
EM oc T a ) between the temperature of peak emi ssion (typically around 3-5 MK) and 1 M K, 



with observed values ranging from 1.70 to 5.17 



2011 



Tripathi et al. 



2011 



Warren et al. 



2012 



(IWarren et al. 



2011 



Schmelz fc Pathak 



Winebarger et al. 



20121 ) . This slope 
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measures the amount of warm plasma (T«l MK) relative to the amount of hot plasma 
(T > 3 MK). Further, the distribution falls off with a very steep slope at temperatures 
above the EM peak; that is, there is very little observable emission hotter than at the peak 



because the emiss ion measure of the p 



ionization effects ( iBradshaw fc Cargill 



asma is low and because of st rong non-equilibrium 



2006 



Reale Orlando 



2008 



In Paper I we found that the large uncertainties in observed a values yield as many as 
77% or as few as zero observed active region EMs that are consistent with low-frequency 
nanoflare heating an d that low-frequency nanoflares cannot explain the upper range of 



observed EM slopes. 



Warren et al. 



(120 111 ) found that very high-frequency heating leads 



to hot, isother mal EMs, which have no t been observed. Warm em ission is also present in 



AR cores (e.g. 



Viall k Klimchuk 



2011 



2012 



Warren et al. 



20121 ). The key to obtaining 



steeper slopes is to enhance the amount of hot emission relative to the amount of warm 
emission, but one must still account for the presence of warm emission. Furthermore, the 
hot component of the EM is itself not generally perfectly isothermal and can extend into 
temperature bins either side of the peak. Increasing the frequency of heating events will 
yield steeper EM slopes, but it cannot be increased so much that the resulting EM is 
effectively isothermal, and the duration of heating cannot be so long that no warm emission 
is ever produced. 

These requirements have led us to explore the possibility of what we term "nanoflare 
trains" , which we define as nanoflares that occur on a single magnetic strand and repeat at 
intermediate frequencies. An intermediate frequency is such that the time interval between 
successive heating events is less than the cooling timescale so that another nanoflare occurs 
on the same strand before the loop cools fully, but is not so high that the heating is 
effectively steady. The heating eventually ceases so that the loop then cools and drains. 
Intermediate-frequency heating can maintain the plasma at high temperature for longer 
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than low-frequency heating, which enhances the amount of hot emission relative to the 
amount of warm emission produced when the plasma finally cools, and steepens the EM 
slope. These nanoflare train properties can therefore satisfy the requirements of yielding 
steeper slopes, ensuring that the hot component of the EM has so me intrinsic width by 



Susino et al. 



2010 



allowing a period of cooling between successive heating events (e.g. 
found that EM peaks in ARs indicate conditions in which the inter-event cadence is shorter 
or of order the plasma cooling time), and accounting for the presence of warm plasma by 
allowing the strands to cool fully following the cessation of the nanoflare train. It is well 
known that warm plasma is overdense compared to hydrostatic equilibrium and cooling 



provi des a natural explanation for this observation (j Warren et al. 



2002 



2003|). 



Spadaro et al. 



We note here that steep EM slopes can also be obtained in a scenario of constant 
heating followed by full cooling and draining when the heating is switched off. The 
longer the heating is applied, the more hot emission relative to warm emission and the 
steeper the EM slop e. This scenario is equivalent to a high-frequency nanoflare train (e.g. 



Warren et al. 



20101 ) . While this scenario is feasible, it is not one that we favor for three 



reasons: (1) it is difficult to imagine a scenario in which the power delivered to the plasma 
remains constant and smooth, even if it were continuous we may expect it to fluctuate 
which could look like individual pulses - the very process of releasing energy alters the 
properties of the magnetic field and the plasma; (2) the hot component of the EM would 
appear much more isothermal than it does; and (3) continual, small-scale bursts of activity 
observed at the limits of current instrument resolution support frequent, short timescale 
events. All feasible theories of coronal heating to-date predict that individual magnetic 
strands are heated in an im pulsive manner. T his includes both magnetic reconnection-like 



processes and wave heating (iKlimchuk 



200J) 
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We will present the results of sixty numerical simulations of nanoflare train heating, 
exploring a wide parameter space, in an effort to explain the properties of observed EM 
distributions and we will show that an observational bias in measured a values may arise 
depending upon the EM reconstruction method employed. In Section [2] we describe the 
numerical aspects of this work and in Section [3] we present and discuss our findings. We 
summarize our results, present our key conclusions and discuss directions for future work in 
Section HI 



2. Numerical modeling 

The current work presents the results of numerical calculations performed to investigate 
the feasibility of coronal loop heating by nanoflare trains. The calculations were carried out 
using the HYDRAD code (Bradshaw & Mason 2003; Bradshaw & Klimchuk 2011; Paper I), 
which solves the one-dimensional equations of hydrodynamics that describe the behavior 
of a two-fluid plasma confined to an isolated magnetic strand. The primary equations 
(c onservation of mass, momentum , and energy) and associated assumptions are described 



in 



Bradshaw fc Klimchuk 



(120111 ) . The particular details of the calculations carried out 
here are summarized in Paper I, with the exception of the implementation of heating by 
nanoflare trains which we will discuss now. 

As in Paper I we assume that the heating arises from the impulsive release of energy 
from the magnetic field and so the total amount of energy available to heat the plasma 
must be limited to the amount of free energy in the field. We preferentially energize the 
electrons and make the assumption that not all of the free energy is released from the field 
during a single heating event, but instead during a series of nanoflares: a nanoflare train. 
The free magnetic energy density is given by 



-7- 



f - {eBp)2 m 
Eb ~ "ST' (1) 

where _B P is the potential component of the field and e parameterizes the level of stress such 
that B s = eB p is the st ress component. We choose e = 0.3 (Dahlburg et al. 2005; Paper I). 



Mandrini et al. 



(|2000| ) studied how the average values of B and B 2 depend on the length 
of the field line, 2L. Table [1] gives results for several observed active regions based on 
Equation 9 of their paper. 

Since -B mm and -Bmax can vary quite substantially we have assumed that B p m B m g 
in our calculation of the free magnetic energy density in the fifth column of Table [TJ This 
is the total energy available for release during the nanoflare train. To determine the peak 
volumetric heating rate of energy release Eh for each nanoflare of the train we need to 
know the length of the field line 2L (so that we can find -Bavg and hence Eb), the number 
of nanoflares N comprising each train, and the duration of each individual heating event 

We have chosen to explore a parameter space in which 2L = [40, 80, 160] Mm, 
iV = [5, 10, 15, 20] and th = [60, 180, 300] s, and carry out 3 sets of 20 numerical 
experiments; 1 set for each value of 2L. The range of t# was decided by considering what 

Table 1: Average magnetic field strength (<S 2 > 1//2 ), free energy and volumetric heating 
rate versus loop length. 



2L 


^min 


-Smax 


-Bavg 


Eb 


Eho 


[Mm] 


[G] 


[G] 


[G] 


[erg cm -3 ] 


[erg cm -3 s 


40 


83 


189 


136 


66.23 


0.03680 


80 


42 


150 


94 


31.64 


0.01760 


160 


18 


89 


51 


9.31 


0.00517 
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Fig. 1. — The temporal heating profile of the nanoflare train for Run 21, where 2L = 80 Mm, 
N = 5, r H = 60 s and t c = 60 s. 
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constitutes an impulsive heating event. The most straightforward way to define impulsive 
heating is to set an upper limit to the timescale such that th « Ac, where Ac is the 
total cooling timescale (thermal, radiative, enthalpy-driven) at the cessation of heating. 
Simulations of the secondary instability of electric current sheets indicate heating timescales 
of order 100 s, though this is highly variable depending on the thickness of the sheet 



(IDahlburg et al. 



20051 ) . For the range of loop lengths, temperatures and densities typically 
encountered in active regions, Ac is in the region of a few thousand seconds (Table EJ). We 
have set what we consider to be a reasonable limit to impulsive heating of about 10% of 
this timescale and investigate a range of r# within this limit. However, we note here that 
longer duration heating may be appropriate for the very longest loops (160 Mm) in our 
study and the consequences of increasing th for longer loops can easily be inferred from 
our results for shorter loops. We chose the range values of N in order to create nanoflare 
trains with total durations Ah that are both shorter and longer than Ac so that we could 
investigate in detail the relationship between the amount of hot plasma (dependent on Ah), 
the amount of cooling plasma (dependent on Ac) and the EM slope. Another factor in the 
choice of N is that the loop lifetime should be consistent with observed lifetimes on the 
order of hours. We have also chosen to keep Eho constant for each set of 20 experiments, 
so that we can focus on the influence of N, th and tq (the time interval between successive 
events) on a (where EM(T) oc T a ). Assuming that the free energy is divided more or less 
equally among each of the individual nanoflares and the temporal profile of the heating is 
triangular (a linear rise and decay, e.g. Figured]), the total energy release per unit volume 
can be found from 



E H = \Nr H E m . (2) 
Setting Eh = Eb and choosing N = 20 and th = 180 s, we find the values of Eho listed 
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in the final column of Table [TJ This choice of Eho ensures that for th = 300 s and all N, 
-^min < B < B max (where B is the field strength implied by E H ). 

The loop geometry is semi-circular with a 2 x 10 4 K chromosphere of depth 10 9 cm 
(several scale heights) at each footpoint, providing a reservoir of mass and energy that will 
be supplied to the corona upon heating. Figure H] shows the initial temperature and density 
profiles for the set of numerical experiments where 2L = 80 Mm, which were found by 
integrating the hydrostatic equations from foot-point to foot-point. The initial conditions 
are energetically negligible compared with the magnitude of energy released during the 
nanoflare train. 

The forward modeling aspect of this work is described in detail by Bradshaw & 
Klimchuk (2011) and in Paper I, though we summarize the salient information here. A 
multi-stranded loop is created such that each strand represents a one second interval of the 
entire heating and cooling cycle. For example, if a heating and cooling c ycle took 3600 s 

loop would consist of 3600 individual strands (e.g. 



hour) then our 


mode 


Guarrasi et al. 


2012) 



Testa et al. 



2005 



procedure provides a reasonable representation of the average state of the loop atmosphere 
when observed. The EM is then calculated in two separate ways. We first calculate an EM 
that has no dependence on instrumental or line-of-sight constraints by directly evaluating 
the integral (EMc = J n 2 dr). This is the true, or model, EM. We can do this because we 
have access to the values of density and temperature for each strand, calculated by our 
numerical model, a luxury not available to observers. We calculate the apex EM (to avoid 
foot-point / moss contamination) and sum over all of the strands to find the total EM. 

We also calculate a synthetic EM along the line-of-sight that Hinode-EIS would see 



if it were to observe loops in an active region that were identi c al to our mode 



use the Pottasch method (IPottasch 



1963 



Jordan et al. 



1987 



Tripathi et al. 



loops . We 



20111) and 
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Fig. 2. — The initial hydrostatic temperature and density profiles for the set of numerical 
experiments where 2L = 80 Mm. 
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Table 2: The 30 spectral lines used to compute the synthetic Hinode-EIS emission measure. 



Ion 


\ A / o T ro 1 an itt ri 

vv aveieng in 


lno- T 
1O gl0 L 


Ion 


\ A / o T ro 1 /ill itT ri 

vv dveieng in 




Mo- V 
IVlg V 


z i o.o t y 


K A K 
0.40 


T?q ytt 
re A.11 


1 Q9 

iyz.oy4 


p. on 
o.zu 


Mo- VT 
IVlg V 1 


9fi8 QQ1 

zoo.yy i 


0.00 


TTo YTT 
re a.11 


iyo. i iy 


9n 
o.zu 


Mo- VT 
IVlg V 1 


97D ^01 

z i u.oy i 


0.00 


Tie. yttt 
re a.111 


ZUZ.U44 


P. 9^ 

o.zo 


q; vtt 

ol V 11 


Z I 0.O04 


O.oU 


TTq "VTTT 

re A.111 


ZUO.oZo 


9^ 

o.zo 


Mg VII 


278.404 


5.80 


Fe XIV 


264.790 


6.30 


Mg VII 


280.745 


5.80 


Fe XIV 


270.522 


6.30 


Fe IX 


188.497 


5.85 


Fe XIV 


274.204 


6.30 


Fe IX 


197.865 


5.85 


Fe XV 


284.163 


6.35 


Si IX 


258.082 


6.05 


S XIII 


256.685 


6.40 


Fe X 


184.357 


6.05 


Fe XVI 


262.976 


6.45 


Fe XI 


180.408 


6.15 


Ca XIV 


193.866 


6.55 


Fe XI 


188.232 


6.15 


CaXV 


200.972 


6.65 


Si X 


258.371 


6.15 


Ca XVI 


208.604 


6.70 


Si X 


261.044 


6.15 


Ca XVII 


192.853 


6.75 


S X 


264.231 


6.15 


Fe XVII 


269.494 


6.75 
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construct an apex EM for each strand (summed to find the total EM) from 30 spectral 
lines chosen from published observational studies. These lines are formed across a wide 
range of temperatures and are listed in Table |5J The forward-modeling procedure by which 
the line intensities w ere synthesized using the conce pt of a virtual detector is described in 



detail in Section 3 of 



Bradshaw fc Klimchuk 



(120111 ) and summarized in Appendix |A] A set 



of EM-loci curves can then be derived by dividing each line intensity by the contribution 



funct ion (calcula ted from a t omic data provided by the Chianti database: 



2009|). Following 



Pottasch 



Dere et al 



1997 



( 119631 ) the values for the EM assigned to each line are estimated 
by assuming that the contribution function is a square function having a constant value 
over a temperature range of width AlogT = 0.3. We have used a density of 10 9 cm -3 in the 
contribution functions because we showed in Paper I that this value brings density sensitive 
lines into better agreement with their neighbors. 

The EM derived in this manner is subject to some of the same constraints as a real 
observed loop and can be compared to published EMs derived from observational studies. 
By calculating these two EMs, we can see what features of the true EM are reliably 
reproduced in the synthetic one and what information about the state of the plasma is 



lost. We can t 
observed EM. 



ren establish a level of confidence in the information that we extract from an 



Testa et al. 



(120121 ) have carried out an extensive study of this important 
issue by comparing known model EMs from a 3D model of a quiet Sun region, with EMs 
derived by tracing different lines-of-sight through the computational domain and treating 
the summed emission as though it were real data. They used the Extreme-ultraviolet 
Imaging Spectrometer (Hinode-EIS) response functions to calculate synthetic spectra and 
the Atmospheric Imaging Assembly (SDO-AIA) response functions to calculate synthetic 
intensities for its wavelength channels. The EMs were then constructed using the Monte 
Carlo Markov Chain (MCMC) method. They found that the EMs derived from EIS 
synthetic data were able to reproduce some characteristics of the model distributions, 
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but showed differences when structures with significantly different density intersected the 



line-of-sight. The EMs derived from AIA synthetic dat a were much 
calculate synthetic EIS spectra in a similar manner as iTesta et al. 



ess accurate. We 



( 12012( 1. by taking 



lines-of-sight through the apex of our target loop. We are primarily interested in AR cores 
and if these can be observed directly from above at disk center then it is reasonable to 
suppose that they are the densest structures along the line-of-sight, and therefore should 
determine the magnitude of the EM. However, we appreciate that this is an important 
observational issue that deserves attention if EMs are to be used as a diagnostic of timescale 
related to coronal heating and we base our estimates of the un certainties assoc i ated w ith 
the E M slope derived from observations on the recent work of 



Guennou et al. 



(2012a 



20131) . 



3. Results 



Table 3:: The results of 60 numerical experiments for 
heating by nanoflare trains. 



Run # 


2L 


N 








A c 


A C 


logio T peak 


a model 


^observed 


1 


40 


5 


60 


60 


540 


2203 


.245 


6.45 


1.07 


(± 0.10) 


0.91 (± 0.20) 


2 


40 


10 


60 


60 


1140 


1949 


.585 


6.55 


1.61 


(± 0.11) 


1.55 (± 0.29) 


3 


40 


15 


60 


60 


1740 


1866 


.932 


6.55 


1.94 


(± 0.28) 


1.93 (± 0.35) 


4 


40 


20 


60 


60 


2340 


1823 


1.28 


6.55 


2.18 


(± 0.40) 


2.24 (± 0.40) 


5 


40 


5 


60 


300 


1500 


4198 


.357 


6.35 


1.92 


(± 0.34) 


1.61 (± 0.25) 


6 


40 


10 


60 


300 


3300 


3844 


.858 


6.35 


2.78 


(± 0.82) 


2.93 (± 0.36) 


7 


40 


15 


60 


300 


5100 


3776 


1.35 


6.35 


3.30 


(± 1.08) 


3.74 (± 0.41) 


8 


40 


20 


60 


300 


6900 


3771 


1.83 


6.35 


3.59 


(± 1.26) 


4.26 (± 0.44) 


9 


40 


5 


180 


180 


1620 


3767 


.430 


6.55 


1.90 (± 0.18) 


1.98 (± 0.35) 
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Table 3:: The results of 60 numerical experiments for 
heating by nanoflare trains. 



m # 


2L 


N 


TH 




Ah 


A c 


A u 
A C 


logio T peak 


^model 


a observed 


10 


40 


10 


180 


180 


3420 


3453 


.990 


6.55 


9 AQ 
z.^i-y 


(4- 


n ao\ 

U.4tZ 1 


9 7fi 
Z. 1 u 


(ZLZ U.^ i ) 


11 


40 


15 


180 


180 


5220 


3429 


1.52 


6.55 


9 8/1 


i ztz 


U.OZ J 


98 
O.Zo 


(4- n KA\ 
IZC U.O^i-J 


12 


40 


20 


180 


180 


7020 


3436 


2.04 


6.55 


O.U ( 


i ztz 


u.oy ) 


O.UO 


^ zt u.oy j 


13 


40 


5 


300 


60 


1740 


3453 


.504 


6.65 


1 QQ 
i .yy 


(4- 

[ ZJZ 


n 1 (\\ 

U. 1UJ 


1 R8 
1 .Uo 


^ztz u.zu j 


14 


40 


10 


300 


60 


3540 


3292 


1.08 


6.65 


9 ^A 


(4- 

[ ZJZ 


U.OO ) 


9 41 

Z.^i-l 


^ztz u.zo ) 


15 


40 


15 


300 


60 


5340 


3270 


1.63 


6.65 


9 87 
Z.o ( 


( _i_ 

1 zc 


n /i i \ 


9 81 
Z.ol 


(4- D OS"! 
(ZC U.Zo J 


16 


40 


20 


300 


60 


7140 


3284 


2.17 


6.65 


o.uu 


(4- 

[ ZJZ 


n aq\ 


O.U ( 


(zc u.oij 


17 


40 


5 


300 


300 


2700 


3530 


.765 


6.55 


9 A"\ 


(4- 

l ZLZ 


U.OU ) 




('-i- n ai\ 
(zc u.'io j 


18 


40 


10 


300 


300 


5700 


3448 


1.65 


6.55 


D8 
O.Uo 


1 ZLZ 


U.OZ J 


o.oy 


(4- n ^/H 

IZC U.041J 


19 


40 


15 


300 


300 


8700 


3446 


2.52 


6.55 


1 AO 


1 ZLZ 


U.UO ) 


41. 10 


(4- n ^ch 
^ zt u.oy ) 


20 


40 


20 


300 


300 


11700 


3457 


3.38 


6.55 


O.UO 


(4- 
{ ZLZ 


n 79"! 

U. 1 Zij 


ZL ^fi 
I.JO 


(4- n «qa 
^ztz u.uo j 


21 


80 


5 


60 


60 


540 


4333 


.125 


6.35 


1 9*? 

i .zo 


I _|_ 
I ZLZ 


U.UU ) 


n 88 

U.OO 


/-i- n 1 9"i 

II U. 1Z J 


22 


80 


10 


60 


60 


1140 


3939 


.289 


6.65 


1 1 1 
1.11 


I ZLZ 


n 1 9"i 

U. 1ZJ 


1 AQ 
i .^ty 


(4- n 1 7"\ 

(ZC U. 1 1 J 


23 


80 


15 


60 


60 


1740 


3699 


.470 


6.65 


1.53 


(± 


0.07) 


1.77 


(± 0.22) 


24 


80 


20 


60 


60 


2340 


3538 


.661 


6.65 


1.77 


(± 


0.14) 


1.97 


(± 0.26) 


25 


80 


5 


60 


300 


1500 


4198 


.357 


6.45 


1.29 


(± 


0.06) 


1.15 


(± 0.13) 


26 


80 


10 


60 


300 


3300 


3844 


.858 


6.45 


2.16 


(± 


0.36) 


2.10 


(± 0.22) 


27 


80 


15 


60 


300 


5100 


3776 


1.35 


6.45 


2.57 


(± 


0.59) 


2.76 


(± 0.30) 


28 


80 


20 


60 


300 


6900 


3771 


1.83 


6.45 


2.84 


(± 


0.73) 


3.26 


(± 0.35) 


29 


80 


5 


180 


180 


1620 


3767 


.430 


6.55 


1.61 


(± 


0.04) 


1.81 


(± 0.22) 


30 


80 


10 


180 


180 


3420 


3453 


.990 


6.65 


2.15 


(± 


0.16) 


2.24 


(± 0.21) 


31 


80 


15 


180 


180 


5220 


3429 


1.52 


6.65 


2.49 


(± 


0.28) 


2.69 


(± 0.26) 
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Table 3:: The results of 60 numerical experiments for 
heating by nanoflare trains. 



an# 


2L 


N 


TH 




Ah 


A c 


A u 
A C 


logio T peak 


^model 


^observed 


32 


80 


20 


180 


180 


7020 


3436 


2.04 


6.65 


9 79 


(4- 


u.oo j 


^ nn 
o.uu 


^ZC U.OU ) 


33 


80 


5 


300 


60 


1740 


3453 


.504 


6.65 


1 P.P. 
1 .00 


i zt 


n r\7\ 

U.U ( ) 


l.oo 


(4- f\ -\Q\ 
II U. lo ) 


34 


80 


10 


300 


60 


3540 


3292 


1.08 


6.65 


9 no 
z.uy 


i zt 


U.ZO ) 


9 

z.uo 


^ Zt U. 10 J 


35 


80 


15 


300 


60 


5340 


3270 


1.63 


6.65 


9 

Z.OU 


i zc 


U.OO ) 


9 An 


( _i_ n 1 7^1 

^ Zt U. 1 ( ) 


36 


80 


20 


300 


60 


7140 


3284 


2.17 


6.65 


9 AA 


i zc 


n ai\ 

U.40 J 


Z.00 


I _|_ n in"! 

ii u. iy j 


37 


80 


5 


300 


300 


2700 


3530 


.765 


6.55 


1 QQ 

i .yy 


i zc 


U. loj 


9 9D 
Z.ZU 


^ZC U.OU J 


38 


80 


10 


300 


300 


5700 


3448 


1.65 


6.65 


9 ^7 

Z.O ( 


(4- 

1 ztz 


n ^i ^ 

u.oi j 


9 89 


^zt u.zo j 


39 


80 


15 


300 


300 


8700 


3446 


2.52 


6.65 


9 88 


(4- 

1 ztz 


n ao\ 


^ 9fi 
o.zo 


^zt u.oo ) 


40 


80 


20 


300 


300 


11700 


3457 


3.38 


6.65 


3 D7 
O.U ( 


( _l_ 

1 zc 


n aq\ 
u.^ty j 


c;7 

0.0 1 


^ZC U.O 1 J 


41 


160 


5 


60 


60 


540 


6612 


.082 


6.35 


n 7Q 
u. i y 


I zc 


U.U 1 J 


l.lo 


^ Zt U.OO J 


42 


160 


10 


60 


60 


1140 


6585 


.173 


6.45 


1 97 


(4- 

[ ZJZ 


n 1 A\ 

U. 14 I 


1 97 


(4- n oi a 
^zt u.zij 


43 


160 


15 


60 


60 


1740 


6577 


.265 


6.65 


1 no 
r.uy 


i zc 


U. 10 ) 


1 £1 

l.Ol 


f-i- n 91 "\ 

I ZC U.Z1J 


44 


160 


20 


60 


60 


2340 


6579 


.356 


6.65 


1 A 1 

1.41 


i zc 


n 1 1 

U. 1 1 J 


1 89 
i.oz 


(4- n 91 "\ 

^ZC U.Z1J 


45 


160 


5 


60 


300 


1500 


6006 


.250 


6.35 


0.96 


(± 


0.07) 


1.22 


(± 0.67) 


46 


160 


10 


60 


300 


3300 


6099 


.541 


6.45 


1.62 


(± 


0.11) 


1.50 


(± 0.23) 


47 


160 


15 


60 


300 


5100 


6112 


.834 


6.45 


2.05 


(± 


0.29) 


2.04 


(± 0.26) 


48 


160 


20 


60 


300 


6900 


6082 


1.13 


6.45 


2.35 


(± 


0.43) 


2.46 


(± 0.29) 


49 


160 


5 


180 


180 


1620 


6718 


.241 


6.55 


1.32 


(± 


0.14) 


1.63 


(± 0.21) 


50 


160 


10 


180 


180 


3420 


6567 


.521 


6.65 


1.78 


(± 


0.06) 


2.11 


(± 0.24) 


51 


160 


15 


180 


180 


5220 


6434 


.811 


6.65 


2.12 


(± 


0.15) 


2.25 


(± 0.19) 


52 


160 


20 


180 


180 


7020 


6425 


1.09 


6.65 


2.35 


(± 


0.21) 


2.55 


(± 0.22) 


53 


160 


5 


300 


60 


1740 


6762 


.257 


6.65 


1.37 


(± 


0.14) 


1.94 


(± 0.21) 
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Table 3:: The results of 60 numerical experiments for 
heating by nanoflare trains. 



Run # 


2L 


N 


th 




A H 


Ac 


A C 


logio T peak 


a model 


^observed 


54 


160 


10 


300 


60 


3540 


6463 


.548 


6.65 


1.86 (± 0.07) 


1.64 


(± 0.13) 


55 


160 


15 


300 


60 


5340 


6344 


.842 


6.65 


2.13 (± 0.15) 


2.02 


(± 0.13) 


56 


160 


20 


300 


60 


7140 


6282 


1.14 


6.65 


2.30 (± 0.24) 


2.28 


(± 0.15) 


57 


160 


5 


300 


300 


2700 


6606 


.409 


6.65 


1.55 (± 0.09) 


2.02 


(± 0.23) 


58 


160 


10 


300 


300 


5700 


6530 


.873 


6.65 


2.22 (± 0.14) 


2.38 


(± 0.20) 


59 


160 


15 


300 


300 


8700 


6471 


1.34 


6.65 


2.56 (± 0.25) 


2.81 


(± 0.25) 


60 


160 


20 


300 


300 


11700 


6490 


1.80 


6.65 


2.79 (± 0.32) 


3.12 


(± 0.28) 



We performed 60 numerical experiments for nanoflare train heating within the 
parameter space described in Section |2j The results of these experiments are summarized 
in Table |3j As an example consider Run 21, for which the temporal profile of the enegy 
release by a nanoflare train is shown in Figure [U The nanoflare train consists of N = 5 
separate heating events of period th = 60 s each, with an interval of tq = 60 s between 
each event. The total heating time Ah is the time from the onset of the first heating event 
until the cessation of the final heating event. The total cooling time Ac is the time from 
the end of the final heating event until the electron temperature falls below 10 5 K. The 
temporal evolution of the average coronal electron and ion temperatures, and the average 
electron density, are shown in the upper two panels of Figure |3j The model and synthetic 
observed EMs derived from the results of Run 21 are shown in the lower panel of Figure El 
where the diamonds are the model values and the plus-signs are the synthetic quantities 
derived using the Pottasch method. EM loci curves are shown as dotted lines. The EM 
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Fig. 3. — Results for Run 21. Upper-left panel: time-evolution of the average electron 
and ion temperatures. Upper-right panel: time-evolution of the average coronal electron 
density. Averages are calculated over the upper 10% of the loop, centered on the apex. 
Lower-panel: model emission measure (diamonds) and the emission measure derived from 
synthetic Hinode-EIS data (plus-signs). Pixel 42 indicates the apex pixel on the model 
detector. 
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peaks at a temperature Tp ea k in the region of 10 6 - 35 K. A linear regression applied to the 
model and synthetic quantities between T^ eg ^ and 10 6 K yields slopes of a moc [ e i = 1-23 
and ^observed = 0.88. The quantity in brackets that accompanies each slope in Table [3] is 
the one-sigma uncertainty estimate of the slope (found with the standard IDL procedure 
LINFIT). It is important to remember that the same atomic physics quantities are used to 
infer the EM from the synthetic line intensities as are used to derive the intensities from 
the model. In the case of actual observations, the adopted atomic physics quantities are 
likely t o be different from the tr ue qu antities, so the uncertainties in the slope are much 



larger ( IGuennou et al. 



2012a 



2013) 



Based on the results presented in Table [3j we can make a number of observations 
concerning the relationship between the parameters that we have explored and the EM 
slopes that we have found. For each set of 20 runs, where 2L remains fixed, it can be 
seen that steeper slopes are obtained with increasing N when th and tq are fixed. In the 
case of Runs [8,20] for iV = 20 we find tt b serve d > 4, which is significantly steeper than 
any of the EM slopes we found in our previous investigation of low-frequency nanoflares. 
This indicates, as expected from the discussion in Section (TJ that sustaining the emission 
close to the temperature of peak emission measure for an extended period is an important 
element of obtaining slopes that are comparable to the upper-range of those that have been 
observed. 

We also find that when considering pairs of Runs where only r# varies (e.g. [1,13], 
[2,14],.. .,[8, 20]), the slope is always steeper for the Run with longer with no exceptions. 
Nonetheless, longer heating timescales for the individual nanoflares of the train are 
not sufficient by themselves to guarantee steeper EM slopes. Consider the group of 
Runs [4,8,12,16,20] for all of which N = 20. Run 4 has the shallowest slope and Run 20 the 
steepest slope, which may fit with expectations. However, Run 8 has a significantly steeper 
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slope than both Runs 12 and 16. It is clear that adjusting m by itself is not sufficient to 
guarantee slopes towards the upper-range of those observed and in consequence we may 
conclude that the slope of the EM is not sufficient by itself to constrain th- 

Further examination of Table |3] yields a connection between the time interval tq 
before the next nanoflare of the train and the EM slope. If we consider pairs of Runs 
for which iV and th are fixed then a pattern emerges. Runs [4,8] and [16,20] both show 
steeper slopes for longer tq- The pattern persists when we include Runs [9-12] in our 
analysis, which have intermediate heating timescales (e.g compare Runs [8,12] and [12,16]). 
Longer time intervals between successive nanoflares consistently yield steeper slopes. This 
relationship between tq and ^observed su gg es ts that while sustaining the emission close to 
the temperature of peak emission measure via larger N and r# is generally necessary for 
obtaining steeper slopes, the interval tq between individual nanoflares is also an important 
parameter when it is shorter than the characteristic cooling timescale of the plasma. This 
may partly be explained by comparing Runs [8,16], [28,36] and [48,56], where each pair of 
Runs has N = 20, th = [60, 300] s and tq = [300, 60] s. The Run with the shorter th and 
longer tq of the pair always has the steeper slope but also a lower log 10 T^ eg ^ and smaller 
A log 10 T between the EM peak and 1 MK. This may have the effect of steepening the 
slope, because the gradient is inversely proportional to A log 10 T. We also expect incr easing 



Susino et al. 



feoid ) who 



tq to effectively smear out the EM in temperature, as found by 
studied the effects of both uniform and localized heating in steady and impulsive regimes. 
They found that extended tq (250 — 2000 s) produced large oscillations in the coronal 
temperature and significant smearing of the EM across a broad temperature range, but we 
confine ourselves to generally shorter timescales for tq (< 300 s) and expect any smearing 
of the EM to primarily affect the hot component. We discuss a method by which tq might 
be diagnosed from this effect toward the end of this Section. 
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Fig. 4. — The ratio A H /Ac versus a moc | e j and ^observed ^ or ^ = ^ ^ ne di amon d- 
signs plot the true (model) emission measure calculations and the plus-signs plot the emission 
measure derived from synthetic Hinode-EIS data. The vertical line indicates the upper-limit 
to a b servec l that can be explained with low-frequency nanoflares. 
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The total duration of the nanoflare train is given by Ah = Nth + (N — 1)tc and the 
cooling time Ac following cessation of the final nanoflare was measured from the numerical 
results. The quantity Ah /Ac, which is the ratio of 'hot' strands (e.g. those affected by 
the nanoflare train) to cooler strands (e.g. those undergoing final cooling and draining), 
is useful for parameterizing the relationship between [N, th, tc\ and ^observed - Fig ure IH 
shows the set of A H /A C plotted against the values of a calculated from the model and 
synthetic emission measures. The vertical line at a = 2.6 in the figure delineates the 
boundary between slopes that are consistent with low-frequency nanoflares (a < 2.6) 
and those that are not (a > 2.6). The general trend that emerges is for larger values of 
Ah /Ac to yield steeper EM slopes. Figure ES also shows that slopes shallower than 2.6 
may be consistent with both low-frequency nanoflares and nanoflare trains with relatively 
low N. The case of low-frequency nanoflares is essentially the limit N = 1. We can now 
understand why increasing tq leads to steeper EM slopes; longer tq increases Ah relative 
to Ac- One further trend that emerges in Figure H] is for an increasing discrepancy between 
"model an d "observed as ^ ne s l°P e steepens, where the slope of the EM derived from 
synthetic Hinode-EIS data increasingly over-estimates the slope of the model EM. Since 
we have followed a forward modeling procedure to calculate ^observed then this suggests 
that the slopes found from EMs derived from real observational data might be subject to 
the same bias that underlies this discrepancy, leading to an over-estimate of slopes at the 
upper-range. 

We now turn to considering the effect of changing 2L on the slopes obtained within our 
chosen nanoflare train parameter space. Comparing triplets of Runs for which N, th and tc 
remain fixed and only 2L varies, e.g. Runs [4,24,44], we find that Observed * s consistently 
shallower for longer loops. Figure [5] shows the Ah/ Ac versus a plots for 2L = [80, 160] Mm 
where it becomes clear that the range of a is more restricted with increasing loop length. 
The steepest slopes for 2L = [80,160] Mm are ^observed = [3-57,3.12]. Figure |5] also 
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shows that an increasing proportion of the slopes yielded by nanoflare trains are consistent 
with the range of slopes yielded by low-frequency nanoflares. This is particularly evident 
for 2L = 160 Mm, where all but two values of observed are below ^-6. We can also see 
that the discrepancy between the slopes of the model EM and the EM from synthetic 
Hinode-EIS data persists for longer loops as the slope steepens. A comparison of the data 
points in Figures 0] and \5\ shows that this property does not appear related to A h /Aq 
because Ah/ Ac ~ 3.5 in Figure @] yields a larger difference between a moc i e i and ^observed 
than A H /A C ~ 3.5 in the upper panel of Figure [5) 

We now discuss the physical reason for the dependence of the EM slope on the 
parameter Ah/ Ac- Returning to Figure [3] for Run 21 (N = 5, r# = 60 s, tq = 60 s) we see 
from the upper panels that the coronal density does not reach its peak until significantly 
after the nanoflare train has ended. Consequently, each strand of the loop is well into its 
final cooling phase at the time of peak EM and the data points from the temperature of 
peak EM to 1 MK are dominated by these cooling strands, which means that there is no 
significant contribution to the EM from the hot strands that are undergoing nanoflare train 
heating. The slope is therefore subject to the analytical limit a ma x determined in Paper I. 

Now, examine Figure |6] for Run 40 (N = 20, r H = 300 s, t c = 300 s), which is at 
the opposite extreme to Run 21 for the set of Runs where 2L = 80 Mm. It can be seen 
from the upper panels that the coronal density now reaches and maintains its peak during 
the nanoflare train, which means that we can expect a significant contribution to the EM 
from the heated strands. The EM shown in the lower panel confirms our expectation. 
The temperature of the hot strands varies between 3.5 MK and 5 MK, with an average of 
4.25 MK, and the peak of the EM is substantially enhanced between between 10 6 ' 6 K and 
10 6 ' 7 K and log 10 (4.25 x 10 6 ) = 6.62. The temperature variation allowed by tq during the 
nanoflare train is Alog 10 T = log 10 (5 x 10 6 ) — log 10 (3.5 x 10 6 ) = 0.15 and, since the EM 
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Fig. 6. — Results for Run 40. 
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is binned in temperature intervals of 0.1 dex, this is sufficient for the neighboring bins of 
the peak to also be populated. In consequence, the hot component of the EM is not quite 
isothermal and it may be possible to diagnose the parameter tq from its width (as we show 
later in this Section). The EM in the region of the peak is therefore strongly enhanced 
relative to the part of the EM that is due primarily to the cooling (post-nanoflare train) 
strands. Compare this with Figure E] in which no significant enhancement is apparent. It is 
straightforward to see that a linear regression applied between the temperature of peak EM 
and 1 MK in Figure [6] will yield a steeper slope than the same analysis carried out for the 
EM shown in Figure |3J 

As long as tq < Ac, the EM of the cool component is set by Ac and the EM of the hot 
component is set by A H , which in turn depends on N, t h and t c (A# = Nt h + (N — 1)tc). 
A longer timescale A# means that a greater proportion of the total number of loop 
strands are at some stage of the nanoflare train. This enhances the part of the EM that is 
determined by the nanoflare train relative to the part that is cooling following the cessation 
of the train. From the perspective of an individual strand, the plasma spends more time 
in the hot state and this translates into a larger EM. We can steepen the EM slope by 
changing N, th or re to increase Ah, but Ac is essentially out of our control because it 
depends primarly on the length of the loop. If Tc > Ac, the situation reverts to that of 
low-frequency nanoflares where the slopes are shallower. 

We can now understand why longer loops tend to yield shallower EM slopes between 
the temperature of peak EM and 1 MK. The cooling time Ac tends to be greater for longer 
loops and for the fixed range of Ah that we consider, the ratio Ah / Ac is generally smaller. 
Following the interpretation described above; a smaller proportion of the total number of 
strands will now be at some stage of the nanoflare train and so the hot part of the EM will 
not be so enhanced relative to the post-train part of the EM. Figure [7] shows the EMs for 
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Fig. 7. — The emission measure distributions for Runs 20 (upper-panel) and 60 (lower-panel). 
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Runs 20 (2L = 40 Mm) and 60 (2L = 160 Mm), where N = 20, r H = 300 s, t c = 300 s in 
both cases. It can be seen that the hot part of the EM is greatly enhanced relative to the 
cool part in Run 20, whereas the enhancement isn't so pronounced in Run 60. In the case 
of Run 60, the most extreme of the 2L = 160 Mm set of Runs, the ratio Ah/ Ac = 1-80 
compared to values over 3 for shorter loops. The enhanced hot part of the EM clearly 
leads to steeper slopes and, therefore, in order to obtain steeper slopes for longer loops the 
quantity Ah must be increased relative to Ac- 

Table 0] shows the distribution of observed EM slopes cool-ward of the peak, 
1.70 < a < 5.17 (Table 3 of Paper I), where the uncertainty in the slope Aa = ±1.0 
(Guennou et al. 2012a,b; Paper I). Given our predicted range of 0.88 < ^observed — 4-56 
and the range of uncertainty, we can see that the nanoflare train heating scenario is 
consistent with 86% to 100% of the observed range of EM slopes. This may be compared 
with 0% to 77% in the low-frequency nanoflare scenario found in Paper I. We encourage 
the gathering of many more observational measurements of EM slopes from large-scale 
surveys of active regions in order to improve these statistical estimates. Furthermore, EM 
reconstructions applied to a far larger set of real observational data, such as would be 
obtained from large-scale surveys, would directly address the question of the ubiquity of 
a hot component to the EM and thereby provide direct evidence either for or against the 
heating scenario described here. 

Given a sufficiently large sample of EM observations it may be possible to estimate 
Tc from the characteristic width of the hot component in the following way. The cooling 
timescales by thermal conduction and radiation are given by 



r cond = 4 x 10 and r rad = ^&=i , ( 3 ) 

where x an d b are given by piece-wise power-law fits to the curve of the radiative loss rate 
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as a function of temperature (e.g. iKlimchuk et al. 1120081 ). and the ratio of these cooling 
timescales is given by 



= ( 9 . 66 x 1Q 5) x ( nL f T b-7/2^ (4) 

r rad 

which can be calculated from an observed EM where T = T^ eg ^ if L and n are known. 
The temperature and density range during the nanoflare train in Figure El togeth e r wit 



Klimchuk et al. 



fl2008h . 



the loop length 2L = 40 Mm and the values of \ an d b given by 
yields r conc l / r rac i = 0.03. We conclude that thermal conduction is the dominant cooling 
mechanism during the interval Tc, in this example, and therefore primarily responsible 
for setting the width of the hot component. Using the temperature evolution in the case 
of non-evaporative conductive cooling; (since the dens i ty do es not change much during 



relatively long trains) given by 



Antiochos fc Sturrock 



19761 ) 



where j3 = —2/5, we find 



T(t) = T 1 



r cond 



(5) 



dT 



I3T 



r cond ^ 

We can approximate the width of the hot component AT by the change in temperature 
during the interval tq- 



(6) 



AT 



/3T 



T c r cond + r ^ 

Since the temperature is usually binned in units of log 10 T we can write 



(7) 
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Alog 10 T^ f T ° (8) 

(In 10) (r cond + r c ) 

The quantity A log 10 T < because the plasma is cooling. Finally, the interval tq can be 
found in terms of the approximate width of the hot component from 

^"^(^^(Alog.oT)" 1 ) ' (9) 

When Tc » T con( j we note from Equation El that Alog 10 T = — 2/5(lnl0), which sets 
an upper limit to the component of the width of the hot part of the EM that can be 
due to cooling by thermal conduction. To diagnose r c in this regime using Equation HI 
temperature bins that are substantially narrower than this limit of |Alog 10 T| « 0.17 dex 
are required and we conclude that the typical bin width of 0.1 dex is too coarse to be useful. 

In the regime where r conc l / r rac i > 1 radiative cooling dominates during the interval 



Tc- U sing the temperature evolution in the case of radiative cooling given by 



(1995) 



Cargill et al. 



1-6 



T(f)=T (l-(l-&) — ) (10) 



r rad- 

and following a similar analysis to that presented for cooling by thermal conduction, we find 



Alog in T^ , — (11) 

(In 10) ((1- b) r c -r rad ) 1 ) 



and 



VT^l-i)- I ) . (12) 



31 



When (1 — b)rc » r rac j we note from Equation [TT1 that Alog 10 T is a positive quantity 
and yet we require it to be negative for cooling. The only way to satisfy this constraint is to 
ensure that (1 — b)r c < T ra( j and this implies r c < r rac i/(l — b), which sets the upper limit 
to the inter-nanoflare timescale tq that can be diagnosed using Equation [12j For example, 
b = —1/2 is commonly used in analytical approxim ations to the radiative loss curve in 



the range log 10 T > 5.0 ( jBradshaw fc Cargill 



20101 ) and this sets the limit t c —> 2r ra( j/3 



as Alog 10 T — > oo. In this regime, we conclude that temperature bins of 0.1 dex may be 
adequate depending upon the relative values of tq and i~ rac [- Following a different approach, 
it may also be possible to constrain the parameter space of N, th and tq from time-series 
st udies of active region emission a cross a range of temperatures, such as those carried out 



by 



Viall fc Klimchuk 



(12011 



20121). 



Figures @] and show that a ODservec i an d a mo( [ e \ diverge as the EM slope steepens, 
such that the slope of the EM derived from synthetic Hinode-EIS data increasingly 
over-estimates the slope of the model EM (although in many cases each slope value falls 
within the error bars of the other one). This is particularly noticeable for a > 3. It may 
be attributed to the estimate of the electron density used in the contribution function 
G(n,T) to derive the loci-curves (EM(T) oc I /G(n,T)). The line intensities in our forward 
model are computed using contribution functions that are density dependent. However, 
when calculating the EM using the Pottasch method the line intensities are divided by 
contribution functions that are no longer density dependent (assuming a fixed value of 
n = 10 9 cm -3 ). The hotter strands have the highest densities, but fixing the contribution 
function density at 10 9 cm -3 for density sensitive lines leads to an underestimate of the true 
contribution function and an overestimate of the EM at higher temperatures. The density 
sensitive lines therefore introduce a bias toward steeper slopes. To mitigate this effect 
we recommend that density measurements accompany EM calculations, which provides 
another good reason for measuring n in addition to needing it for estimates of Ac. 
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The shape of G(n,T) can also play a role in the discrepancy between the model and 
predicted slopes because the Pottasch method assumes that G(n, T) is constant in the range 
A log 10 T = ±0.15 to either side of the peak formation temperature of the emission line, and 
zero everywhere else. The constant value is equal to the average of the actual G(n, T) in 
this range. If two emission lines with the same peak formation temperature have G(n, T)'s 
with appreciably different shape (e.g. one is substantially narrower than the other) then 
this could also affect the more density sensitive lines, leading to a discrepancy between the 
model and predicted values. 



4. Summary and Conclusions 

We have run a series of numerical experiments to explore coronal loop heating by 
nanoflare trains with the intention of extending our previous study, presented in Paper I, of 
heating by low-frequency nanoflares. The aim of the present study is to determine whether 
nanoflare trains are able to overcome the intrinsic limit to the EM sl opes (from the EM 



peak to 1 MK) that we demonstrated in our previous work (see also: 



Mulu-Moore et al. 



2011) so that observationally measured slopes of a > 2.6 can be explained. By constructing 
EMs from synthetic Hinode-EIS data, using spectral lines formed over a wide range of 
temperatures, we predict slopes in the range 0.88 < a ODservec l < 4.56 for the nanoflare 
train heating scenario encapsulated by the extent of the parameter space that we have 
chosen for our study. 

Though we have not attempted to match any specific observations in the work 
presented here, we have reached a set of conclusions concerning the broad properties that 
the nanoflare trains must possess to be consistent with observations. A number of our Runs 
yield temperatures of the peak EM that are probably too low to agree with the observed 
range in active region cores; for example, Runs [5-8,21,41,45] have l°gio^peak = 6-35. 
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We may therefore exclude the particular combination of parameter values that led to this 
result for our fixed value of Eho, which are in general: small N; short r#; and long tq- 
Stronger heating would boost log 10 T^ ea ^ into the required range, but must of course 
remain consistent with the free magnetic energy. 

In the case of loops of fixed length 2L we find steeper EM slopes with increasing N, 
th and re- For fixed Ejjo the parameters N and th determine the total energy released 
during the nanoflare train and tq determines the width of the hot component of the EM. 
We have parameterized the relationship between [N, th,tc ] and ^observed usm g the single 
parameter An/ Ac, which is the ratio of the number of strands in some stage of nanoflare 
train heating to the number of strands undergoing cooling and draining following the 
cessation of the nanoflare train (Figures H] and [5]). We found that the general trend is for 
^observed ^° mcrease with Ah/ Ac- The physical reason for this is that the EM of the cool 
component is set by and the EM of the hot component is set by An = Nth + (N — 1)tc- 
A longer timescale Ah means that a greater proportion of the total number of loop strands 
are in some stage of heating. This enhances the part of the EM that is determined by the 
heating relative to the part that is cooling and draining. 

In the case where 2L is allowed to vary we find it relatively easy to obtain slopes 
consistent with the upper-range of those observed for shorter loops (e.g. 2L = 40 Mm, 
Figure HJ), within the parameter space of nanoflare train properties that we explored. The 
range of slopes for longer loops (2L > 80 Mm) that we found within this parameter space 
are largely consistent with low-frequency nanoflares. All but two values of a ir)g ervec i 
are consistent with low-frequency nanoflares in the case of 2L = 160 Mm. However, 
we also know from our investigations that obtaining steeper slopes for longer loops is a 
matter of allowing the ratio Ah/ Ac, and therefore the quantity Ah, to increase. Given 
an observationally measured EM slope of ~ 3, Figures H] and [5] show that the expected 
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ratio Ah/ Ac ~ 1-5 — 2.0 but we cannot determine Ah uniquely from this information 
alone. However, given the length and density (from spectral line ratios) of the structure 
then it is possible to estimate A q as the radiative and enthalpy-driven cooling timescale 



(e.g. 



Bradshaw fc Cargill 



20101 ) from the temperature of peak emission measure (a known 
quantity from the EM plots). Alternatively, it may be possible to measure Ac by following 
the cooling structure via emission from recombining ions or successively cooler wavelength 
channels. We urge observers to include such estimates (2L, n and / or Ac) along with all 
measurements of EM slopes, in order that possible heating scenarios may be more rigorously 
tested against observational data. 

We have found that the properties of 86% to 100% of active region cores can be 

consistent with heating in the form of nanoflare trains, where we have assumed that not all 

of the free magnetic energy is released at once, but rather in a series of individual nanoflares 

that we refer to as a train. However, there are alternative scenarios. In one example, all of 

the energy available for heating may be released and the magnetic stress then builds up 

again in the interval between nanoflares. In a modeling sense this is similar to the scenario 

that we have explored here, but with a smaller amount of free energy in the field before 

the onset of heating because it only has to supply a single event before recharging. The 

parameter tq then becomes a recharging timescale for the magnetic field. Suppose that 

the storage of energy in the field is a consequence of the field lines becoming twisted and 

tangled due to convective motions on the surface and the energy is released when some 

critical condit ion is satisfied (e.g. the angle betwe en neighboring field lines exceeds some 

' ^1 1 



critical value: 



Parker 



1988 : 



Dahlburg et al. 1 120051 ). If the surface motions do not vary too 



strongly then it is reasonable to suppose that the energy stored in the magnetic field, and 
ultimately released, be proportional to the recharging timescale. We cannot predict when 
reconnection will occur, but we can place constraints on this timescale. It cannot be so 
short that the heating is effectively steady, because this would give rise to a more or less 
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isothermal emission measure. Conversely, for EM slopes greater than about 2.6 it cannot 
exceed a cooling timescale because this would be the limit of low-frequency nanoflare 
heating in which slopes are constrained to a range of values less than 2.6. The recharging 
timescale must be of just the right length such that the width of the hot component of the 
EM is consistent with what is observed. 

Previous studies of the relationship between the EM and the timescales related to 
the heating mechanism have led to findings that are co nsistent with o urs, t hough for 
conceptually somewhat different scenarios. For example, 



Testa et al. 



(120051 ) explored a 



parameter space of heating pulses with timescales that are longer (300 — 1100 s) than those 
we have chosen, though the total duration of heating is similar, and with the energy release 
localized at the loop foot-points. When the spatial scale of heating is sufficiently small a 
thermal instability arises, leading to an ongoing cycle of heating and cooling which gives 
rise to the warm component of the EM in their study. Their cooler loops reach similar 
temperatures (T = 10 6,48 K) to some of ours, but their hotter loops reach significantly 
higher temperatures (T = 10 7 K) due to stronger heating, and they adopted a loop length 
of 2L = 200 Mm that is some 40 Mm longer than our upper limit. They find EM slopes 
approaching 5.0, which are steeper than the slopes we find for our 2L = 160 Mm loops but, 



as we have already discussed, extending Ah relative to Ac (consistent with estimat e s of t 



lie 



Testa et al. 



fl2005f ) 



free magnetic energy) will result in steeper EM slopes. In the case of 
their EM slopes are determined by the relative lengths of the heating and cooling phases of 
the cycle that are ultimately governed by the timescales associated with the evolution of 
the thermal instability. 

It is abundantly clear that different sets of observational data to complement EM 
slope and related measurements are required to help constrain the timescales on which the 
active region heating mechanism operates, and ideally to help constrain the mechanism 



-36 - 



Ugarte-Urra fc Warren 



(120121 ). which 



itself. One such data set may include the findings of 
indicate that active region variability evolves over time as do the relative contributions 
to the EM of hot and warm plasma. This suggests that the frequency of heating may 
also change, with ol der ARs heated by h igher-frequency events, and provides a likely 



explanation for why 



Warren et al. 



(120111 ) found that high-frequency heating was most 
consistent with their observations, whereas others have found that low-frequency heating 
is most appropriate in other cases (e.g. Paper I). An alternative data set may be found 
in the temperature dependent pattern of red- and blue-shifts that are observed in active 
regions. An extension to the work presented here will be an analysis of the Doppler-shifts 
as a function of temperature associated with low-frequency and nanoflare train heating in 
order to determine whether they are consistent with the patterns of flows that are observed 
in active regions, and to investigate their potential to provide additional constraints on the 
parameter space of heating properties. This will be the subject of the next paper in this 
series. 
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A. Spectral Line Synthesis 



The calculation of the total (wavelength inte grated) intensity of a lin e in in strument 



units (DN pixel s ) is described in detail by 



Bradshaw fe Klimchuk 



(1201 II ) and 



summarized for convenience here: 



/ = IR(A) x G(X,n,T)x < EM >, 



(Al) 



where IR(A) is the instrument response function (units: DN pixel -1 photon -1 sr cm 2 ), 
which is the product of the effective area, plate scale, and gain of the instrument. The 
Hinode-EIS response functions are available in the SolarSoft package. 

G(X, n, T) is the contribution function, defined as: 



G(X,n,T) 



0.83 x Ab(Y) xYiX e(X,n,T) 



(A2) 



4tt x (hc/X) 

where 0.83 is the protomelectron ratio, Ab(Y) is the abundance (relative to hydrogen) 
of element Y, Yi is the population fraction of charge state % of element Y, e(X,n,T) is the 
emissivit y of the line (units: erg s -1 cm 3 ) calculated using the EMISS_CALC function 



(Chianti: 



Dere et al 



1997 



20091 ). and hc/X is the photon energy (units: erg). 



< EM > is the spatially averaged column emission measure in the pixel (units: cm 5 ): 



EM mx 

< EM >= — P^E, (A3) 
A pix 

where EM^ X = n 2 e <$/ is the emission measure in the pixel and ^4pj x is the pixel area. It is 
straightforward to see how the emission measure can be calculated from the line intensity, 
given the contribution function and the instrument response. 

The intensity per unit wavelength of the spectral line is given by convolving the total 
intensity with the line broadening function: 



„p I I, (A4) 

V 7T<T \ CT 
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where Ao is the rest wavelength of the line and AA is the Doppler shift due to the 
line-of-sight component of the bulk flow, given by AA = ^v. The total line width, er, is 
determined by the thermal width and the instrument width: 



^ = -^ h + 0.36^ ns , 



(A5) 



where 0+^ = ^"u^h anc ^ suitable values for <jj ns are 60 mA (EIS short wavlength channel) 
and 67 mA (EIS long wavelength channel). In calculating the spectral line profile one must 
be careful to ensure that I = f I\d\. 

The line intensity from each cell of the numerical grid are then added to the total for 
the appropriate detector pixel(s). The correct pixel is determined by projecting each grid 
cell onto a single row of detector pixels and establishi ng what proportion of the em ission 



Bradshaw fc Klimchuk 



201 lh . This 



from the grid cell falls onto each pixel (see Figure 1 of 
procedure is repeated for all of the spectral lines of interest, for all of the grid cells of the 
strand, and for every strand comprising the loop. The result is then a prediction of the total 
emission along the loop, due to the contributions from the many sub-resolution strands at 
different stages of evolution, as would be measured by a particular observing instrument. 
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Table 4: Distribution of observed emission measure slopes cool-ward of the peak. 





«-A« 


a 


a + Act 


a < 1.0 


3 








1.0 < a < 1.5 


5 








1.5 < a < 2.0 


3 


3 





2.0 < a < 2.5 


6 


5 





2.5 < a < 3.0 


2 


3 


3 


3.0 < a < 3.5 


2 


6 


5 


3.5 < a < 4.0 


1 


2 


3 


4.0 < a < 4.5 





2 


6 


4.5 < a < 5.0 





1 


2 


a > 5.0 








3 




100% 


100% 


86% 
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